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SUMMARY 


The purpose of this study was to investigate an impedance 
approach for determining low-frequency environments. 

The study consisted of four major tasks: 

1. development of the equations of motion used in the 
impedance technique; 

2. demonstration of the technique using flight data from 
Titan launch vehicles; 

3. application of the technique to an existing Shuttle 
payload ; 

4. development of the criteria and philosophy of the use 
of the technique for future payloads. 

Results indicate that the method investigated is very 
promising. It gives good results in the frequency domain as an 
initial output and transformed time domain results are also good. 
Its use in certain phases of the payload development cycle could 
substantially reduce both the time and cost required in current 
payload loads integration cycles. 



INTRODUCTION 


Analysis of various payload configurations is a very lengthy 
and expensive task. Extensive individual models of the payload 
and booster must be mathematically coupled for the final system 
models and then elaborate time domain response analyses conducted. 
The task of the analysis integration for all the organizations 
involved as well as the manpower and schedule devoted to model 
coupling and loads computations drive these costs skyward. In 
reusable launch vehicles, payloads are expected to exhibit 
numerous variations in configurations, combinations, and experi- 
ments, and the detailed approach of the past will not be a viable 
cost-effective technique. 

A preferable technique would eliminate the necessity for 
detailed coupled models as well as the need for an integra- 
tion task. If possible, the technique would allow the payload 
organization of designers and dynamicists to generate, with 
model information from the booster organization, their own 
payload response and loads predictions. 

This study was conducted to determine the feasibility of 
a new impedance technique for determining payload low fre- 
quency environments. By accounting for the dynamic coupling 
of the payload and booster in the equation of motion in the 
frequency domain, the analytical effort is diminished by elimina- 
ting the final eigensolut ions as well as reducing the equations 
to simple complex transfer function multiplications. In addi- 
tion, the model requirements of the booster consist of free- 
free unloaded interface modal characteristics. Therefore, the 
task of integrating the loads analysis can be accomplished by 



obtaining a set of "standard" booster model data and the payload 
organizations computing their own loads analysis cycles. 

The information presented in this report includes results 
of the use of the impedance technique on Titan flight data as 
well as predictions of the low frequency environments for a 
proposed Shuttle payload. The requirements for implementing 
the impedance technique and it's feasibility are discussed. 


SYMBOLS 


[c] discrete damping matrix 

DFI development flight instrumentation 

FFT fast Fourier transform 

{F„} external force/torque vectors 

E 

{f} interface force/torque vectors 

[i] unity matrix 

[IMP] complex impedance matrix 

j square root of negative one 

[k] discrete stiffness matrix 

LDEF Long Duration Exposure Facility 

LRC Langley Research Center 

[m] discrete mass matrix 

[PADM] complex point admittance matrix 

[cj)] matrix of modal vectors normalized to unity 

generalized mass 

S Laplace variable 

SRB Solid Rocket Booster 

STS Shuttle Transportation System 

SP Spacelab single pallet 

[T] discrete static reduction matrix 
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[tadm] 

A 

OJ 

C 

a 

{q>, (q>, (q) 

iO, (U, W 

{'rurur) 
©, m, io 


complex transfer admittance matrix 

input frequency to modal frequency ratio 

modal radial frequencies 

damping coefficient 

input radial frequencies 

discrete motion of a flexible body 

modal degree of freedom motions of a booster 
model due to external forces 

modal degree of freedom motions of a booster 
model due to interface forces 

modal degree of freedom motion for a constrained 
model 


Subscripts : 


B 

booster 





E 

external 





I 

interface 





IB 

booster mode shapes 

at 

the 

interface 


i 

frequency increment 





FB 

booster mode shapes 

at 

the 

external force 

points 

fb 

booster mode shapes 

at 

the 

interface force 

points 

P 

payload 





1 

payload number one 





2 

payload number two 
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IMPEDANCE TECHNIQUE 

Development of the Equations of Motion 

The current direction in analytical technique improvements 
include efforts in developing low-cost, low frequency (100 Hz 
and below) payload environment prediction techniques. These 
improvements are aimed at taking advantage of the reusable or 
repeated boost vehicle concept. It is hoped that better pro- 
ductions of the expected payload environments based on previous 
flight data will eliminate a major portion of the conservatism 
in payload design and test. Most design or test specifications 
have unnecessary margin due to uncertainties in these environ- 
ments. The cost of generating extensive response analyses 
necessary to define the environments fully are becoming pro- 
hibitive as well as still containing a few areas of concern 
with respect to true vehicle excitation. Not all future payloads 
will be extremely concerned with extensive loads analyses. 
Depending on the weight margins involved, large margins of 
safety could be used in some payload designs to insure payload 
survivability. Primary and/or secondary structure in these 
"light-weight" payloads can be stiffened to force the payload 
modes above the areas of maximum gain in the booster transfer 
characteristics. Areas of concern in the coupled response to 
vehicle dynamics can be analyzed with large margins of safety 
in design to cover uncertainties. 

It is expected, however, that there will still be payloads 
that can be classified as weight critical. These payloads must 
be designed with a goal of optimizing the size of structural 
members with a desired result of demonstrating with analysis 
and/or test the smallest acceptable design margin. Payload 
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loads analysis of this category generally go through approxi- 
mately three phases of analysis. Preliminary loads are cal- 
culated for initial sizing of members. The next two phases 
are iterative on the design since the analysis lags the detailed 
member design. Considering the magnitude of analyses necessary 
for weight critical payloads as well as the number of launch 
and/or landing events required to be considered, it is evident 
that any technique shown to be able to reduce analysis cost 
could have a major impact on design costs. 

Standard techniques used for integrated loads analyses 
include tasks of compiling and coupling the math models from 
all the organizations having primary structural responsibilities 
For example, in the case of a Titan III-E launch analysis, this 
coupling requires data from at least three organizations: two 

boost vehicles and a payload. A majority of the cost associated 
with this procedure lies in two categories. Large expenditures 
are necessary to: 

1. coordinate information necessary for the analysis; and 

2. create the coupled models and perform the required load 
predictions . 

The magnitude of these efforts is probably most evident when 
late program design changes are necessary. Too often at this 
point projects have expended analytical apportioned budgets and 
it becomes either a guessing game of what impact the design 
change has, or the project is forced into a late program costly 
test series. 

The analytical approach studied in this contract is called 
an impedance technique because it deals with frequency domain 
analysis. The objectives considered when developing the tech- 
nique were based on three major points: 


1 . 


a desire to take advantage of standardization of boost 
vehicles in low frequency environment prediction; 

2. a preference for spectral analysis to gain insight 
into the frequency content of the environment; and 

3. a requirement for cost reduction in the areas of 
analysis and model integration using a comprehensive 
payload user's guide. 

These original objectives were based on the assumption that 
flight data from previous flights, sufficiently describing the 
low frequency environments during those flights , would be 
available and could be used to predict the low frequency 
environments for the following flights. This method relies 
on the ability to extract the feedback characteristics of the 
previous payload and add the effects of the payload of interest, 
with the assumption that the booster transfer characteristics 
don't change. With previous flight data and the developed 
impedance technique, payload designers can eliminate cost by 
not only generating low cost environment predictions but save 
additional effort by eliminating an enormous integration effort. 
The desired result of this approach is to rewrite the equations 
of motion of a physically coupled structural system in such a 
manner as to be able to isolate the dynamic feedback character- 
istics of either physical body. When analyzing payload/booster 
systems, the task involves resolving the coupled system into 
the equivalent payload and booster. Once the equations have 
been resolved in this manner, the dynamic feedback of one 
payload can be easily replaced by another without solving for 
a new set of coupled eigensolutions . 

The key to the approach can be found in the description of 
the interface forces between the payload and booster. (In the 
following development, the forces and accelerations referred to 



are shown schematically in Figure 1.) The first step is to 
develop the equations of motion of the booster alone. Consider 
the interface accelerations for this free-free body {q^}, as 
being the sum of the accelerations due to the external forces, 
{Fg}, and those due to the interface forces, {f}. 

Or, 

{q T } = {q T > + {q T ) (i) 

1 1 e 1 f 

By substituting appropriate expressions for the accelerations 
on the right-hand side of (1) , the total interface accelerations 
will account for both the external forces being transferred 
through the booster and the dynamic feedback of the payload. 

Evaluating the term for the accelerations due to the 
external forces first, consider the following general equations 
of motion for the booster with no payload (free-free) . 

[>%]{q B } + [C B ]{q B ) + [Kg]{q B } = {Fg} (2) 

This set of discrete coupled equations for the booster can be 
simplified by modal substitution to yield a set of uncoupled 
differential equations, or 

[l]t5 B J + r2C B to B Jf4 B } + - 0 fb ] T (F e } (3) 

To solve the differential equations in (3) in the frequency 
domain, first assume zero initial conditions and take the 
Laplace Transform of both sides: 

(Cl]S 2 + fcCgWgJS + !>g]){£ B (s)> = [0 fb ] T (F e (s)} (4) 

By taking advantage of diagonal property of the generalized 
mass, damping and stiffness, equation (4) reduces to 

Es 2 + 2 C b w b S + W 2 J{C b (s)} = [4> fb ] T EF e (s)} (5) 

Equation (5) would serve as a start for a displacement trans- 
formation, however, one further substitution yields acceleration 
as a function of the Laplace variable, S. 




Figure 1. Schematic of Typical Payload/Booster System 
Showing Pertinent Accelerations and Loads. 
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( 6 ) 


S + 2C b w b S + w b 


[ ^]{C B (s)> = [^ fb ] T {F E (s)} 


Or, 


{? B (s)> = |~ 


S + 2? B WgS + 


2^^FB-^ { V 8 > } 


(7) 


To convert to the frequency domain, now substitute S = jQ, 




1 


■J[^f B ] T {f e } _ 


( 8 ) 


(^5 - nj) + j(2c B ai B fl i ) 

i = 1,2,3,... no. of inputs frequencies 
Again substituting, 

X i = V W B 

i = 1,2,3,... no. of input frequencies 
Equation (8) becomes 


(9) 




l 


E } ± 


( 10 ) 


(1 - + j(2C B A i ) 

i = 1,2,3,... no. of input frequencies 
Substituting back to the discrete coordinates expressed in 
equation (2) , the final expression for the interface accelera- 
tions, as a function of the input frequency, will be 


~ x l 

{q T } E = [4 > ib ]E— n 

i (1 - a|) + j(2C B A.) 




(ID 


Or, 


i = 1,2,3,... no. of input frequencies 


{ Ve. ■ [tadm].^}. 


( 12 ) 


i = 1,2,3,... no. of input frequencies 
The coefficients matrix, [TADM]^, represents the complex 
transfer admittance from the point of application of the 
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external forces, e.g. engine thrust time histories, to the inter- 
face degrees of freedom accelerations. 

Next, to derive the expression for the interface accelera- 
tions due to the interface forces, {f}, a set of differential 
equations similar to equation (3) can be written as 

[i]{€b> + E2 ?B 0) B ]{|'} + [oj 2]{5 b } = -[<J) fB ]T{£} (13) 

Here the distinction in the booster modal responses is to 
separate those from the external forces (equation 10) , and 
those from the interface forces. Following the previous pro- 
cedure for solving and simplifying the equations of motion, the 
modal accelerations as a function of input frequency will be 


= i- 


x 


’B 


(1 - xp + j(2C B X ± ) 




(14) 


i = 1,2,3,... no. of input frequencies 
Then the discrete accelerations at the interface are written as 


^IB^~ 


x 


f l “ (1 - xp + j(2C B A ± ) 




(15) 


Or, 


i = 1,2,3,... no. of input frequencies 
{qj.} = [PADM]^{f K 


(16) 


i = 1,2,3,... no. of input frequencies 
As in equation (12), [PADM]^, is the matrix of coefficients for 
the complex point admittance for the booster at the interface 
and is the expression for the acceleration response at the 
interface due to the interface forces. 

Substituting equation (12) and (16) into (1) yields 

{q T } . = [TADM].{F> + [PADM] . (f } . (17) 

I x xE..- x x 


i = 1,2,3,... no. of input frequencies 



The effect of the payload on these interface accelerations comes 
through from interface forces, {fK, which in most cases are 
the reactions to the inertial loads of the payload. With this 
in mind, the final step in the derivation involves expressing 
these interface forces in terms of the accelerations of the 
payload (Figure 1). To accomplish this, consider the equations 
of motion of the payload (free-free) . 

[MpHqp} + [C p ]{q p } + [K p ]{q p } = {f } (18) 

In equation (18) , the degrees of freedom for the payload 
include, again, the interface attach point degrees of freedom 
as well as those for points throughout the structure. The 
following substitution can be made to write the motions of the 
payload in terms of the sum of those relative to the interface 
and the interface motions themselves. 

{q p } = {? p } + [T]^} (19) 

If the substitution is made for modal coordinates then (19) 
becomes 

{q p } - C* p ]{? p } + [t] } (20) 

Limiting for the moment, the constraints of the payload model at 
the interface to those resulting in a statically determinate 
interface, the equations of motion from equation (18) can now 
be shown to be 



f 


( 21 ) 



Separating the upper set of equations in (21) , 

[X]{fp} + [2 ?p( o p ]{i p } + = -[<l> p ] T [M p ][T]{q I } (22) 

from the lower set, 

[T] T [M p ][<() p ]{f p } + ([M p ] + [T T M p T]){q I } = {£} (23) 

It is obvious that a substitution for the interface forces, 

{f} , in equation (17) can be obtained. 

Taking the Laplace transformation of equation (22) and 

simplifying as before, we have 

x 2 

(L> = f 5 — J[<f'p] T [M p ][T]{q } (24) 

P i (1-A; . ) + j (2C p X pi ) i 

i = 1 , 2 , 3 ,... no. of input frequencies 
Now substituting (24) into (23) 

X 2 

([T] T [Mj[4> p ]| 2 — C^p] T [Mp]CT: + [Mj] 

L( 1 ' A Pi> + ^ 2 ^p^Pi^J 

+ [T] T [M p ][T]){q‘ I } i = (f) i (25) 

i = 1 , 2 , 3 ,... no. of input frequencies 
or 

[ IMP ] . { q T } = (f} (26) 

1 1 i i 

i = 1 , 2 , 3 ,... no. of input frequencies 

Here, the complex coefficients matrix [iMP]^ represents 
the impedance or "complex inertia" of the payload at the 
payload/booster interface. 

Now all of the complex transfer function characteristics 
have been derived that are necessary to fully describe the 
vehicle coupled response and the final form of the coupled 
impedance can be generated. Substituting (26) into (17) yields 
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I 


{q T } = [TADM] {F } + [PADM] [lMP].{q_} (27) 

-*-i i ^ i i i-*-i 

i = 1,2,3,... no. of input frequencies 
Or, rearranging 

([I] - [PADM] [IMP] ){q T } = [TADM] {F^} (28) 

i i 1 i i ^ i 

i = 1,2,3,... no. of input frequencies 

This is now the final form for the Impedance equation where 
the matrix of complex coefficients for the interface accelera- 
tions represents the coupled impedance of the payload/booster 
combination, and the right hand side of (28) represents the 
complex "psuedo" generalized force. Since the external forces, 

{F >, are in complex spectral form, and the coefficient matrices 
E 

are complex as well, the desired spectral interface accelera- 
tion can be obtained with equation (28) using nothing more than 
complex multiplication. This simplicity greatly reduces the 
cost of generating the low frequency environmental predictions. 

If we now consider some new payload for the same booster 
with the exact same set of external force/torque vectors 
(Figure 2) , then the interface accelerations for the new pay- 
load will be 

(q T ) = ([i] - [padm] [imp] )~ 1 ([i] 

1 2i i 2i 

-[PADM] [IMP] ) (q } (29) 

i li 1 li 


i = 1,2,3,... no. of input frequencies 
given the same arrangement of interface tie-down points. 

The advantages of analyzing the interface environments 


in this manner appear to be significant. The resultant set of 
coefficients for {q^.} in equation (29) act as a coupled 
analytical filter that modifies the environment for payload 


number 1 to arrive at an environment for payload number 2. 
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Payload #2 





O 


Payload 

#2 


Booster 





Mi {f 


Figure 2. Schematic of Impedance Technique for Replacing 
One Payload Feedback with Another. 
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A payload organization could obtain the complex impedance from 
the proposed data bank that corresponds to the previous inter- 
face accelerations and generate new flight predictions for the 
desired interface. Another possibility could be to synthesize 
the necessary transfer characteristics from flight data rather 
than relying solely on analytical model representations. A 
third possibility is to use equation (29) to evaluate design 
changes on a particular payload at a particular set of attach- 
ment locations. Thus, it appears that either equation (28) or 
(29) could be utilized to significantly reduce the life-cycle 
cost of loads analyses. 

At this point, a brief review of the limiting assumptions 
that were used to derive the final equations (28) and (29) is 
important. In generating the payload impedance characteristics 
the assumption used was a statically determinate interface. 

This was done to simplify the complexity of describing the 
force distributions at the interface. This need not be done, 
however, for a set of equations similar to (21) can be derived 
that are somewhat more complex and cumbersome. The velocities 
and displacements at the interface degrees of freedom must be 
tracked in this case in order to solve for the interface 
forces. However, the majority of interface conf igurations are 
statically determinate by design, therefore, this study con- 
sidered only this case. 

The assumptions with the most impact on the technique 
affect the use of equation (29) more than equation (28) . 

Neither are the forcing functions applied to the vehicle 
exactly repeatable nor are the transfer characteristics the 
same for all booster configurations. In the case of a booster 
system like the Shuttle Transportation System (STS) , these 
characteristics change for each location in the Orbiter bay 
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and for different combinations of payloads in multiple payload 
flights . 

Therefore, the effort during this study was directed 
primarily at determining the feasibility of a practical use of 
the technique. Discovery of pitfalls in flight data analysis 
and assessment of the impact of the assumption made were 
necessary goals. 

To assess the accuracy of the technique, specifically the 
calculated results of equation (28), time domain response result 
were compared with those from the impedance technique for a 
simple system (Figure 3) . This evaluation proved beneficial in 
a number of areas. 

The degree of freedom describing the motion of mass M 2 is 
considered the interface motion for this check case, with the 

external forces, F^, applied at M„ . 

E 3 

One discrepancy became apparent from the beginning per- 
taining to the modal damping. Various techniques have been 
tried to obtain the correct damping for a modally coupled 
model. For example, one technique tried on the Viking Program 
involved calculating the triple matrix product of the coupled 
set of modes and the uncoupled diagonal damping matrices and 
discarding the resulting off-diagonal terms to maintain an 
uncoupled set of equations. This assumption does force the 
coupled damping to be a function of the coupled modes but 
throws away, in essence, the coupling effect from mode to mode. 
The problem lies in the manner in which the damping is measured 
and assumed for a payload and booster. If the damping is 
measured for a typical payload in a modal survey, the recorded 
damping is a function of the boundary conditions, for these set 
the relative motions for each degree of freedom. When the 
payload is physically attached to the booster, a new set of 





Figure 3. Schematic of Three Mass Check Model 





boundary conditions exist and, therefore, a new set of relative 
motions or mode shapes also exist. For small damping, less 
than 1%, the previous philosophies about how to generate 
coupled damping may be adequate for determining payload loads 
and environments, but the impedance technique, equations (28) 
and (29), just exemplifies the discrepancy. Since the technique, 
when accounting for the payload feedback in the equations for 
the interface acceleration, is merely replacing the need for 
the coupled eigensolutions , the modal damping used in generating 
the complex impedance of the payload must be one that will 
result in the expected coupled analyses. If the impedances are 
generated assuming 1% damping for all payload and booster modes, 
for example, and the equivalent coupled model has coupled the 
damping and discarded all off-diagonal terms, there will be 
differences in the two transfer functions since some of the 
resulting diagonals of the coupled damping will be larger than 
1%. It should be noted that these large discrepancies due to 
the damping differences may be evident only in the frequency 
domain. In the time domain these differences may make only 
minute errors in the response. 

In order to get comparable transfer functions for both the 
impedance check and the coupled time domain check, the discrete 
dampings were assumed as shown in Figure 3, thus the transfer 
functions were the same regardless. This point is important 
in light of future comparisons of existing loads analysis tech- 
niques and those similar to that developed above. There are, 
however, possible solutions (Reference 1) and the point will 
not be labored further here. 

Use of the Fast Fourier Transform (FFT) in obtaining the 
spectral data also presents some problems as discovered by 
experimenting with the sample problem. The FFT is more commonly 



used for spectral analysis, e.g., auto-spectral densities, 
coherence, etc., of random data (References 2 and 3). Certain 
errors are introduced by signal truncation of time histories. 

In random data analysis, these errors are usually minimized by 
taking a number of averages or statistical degrees of freedom 
(Reference 3) . The effects of truncation on transient or 
complex signals cannot be averaged out since the signals are 
of finite length and do not repeat themselves. These truncation 
errors can be significant in generating Fourier spectra, 
expecially when these spectra are used to drive a system. The 
system transfer function has a certain value at a certain fre- 
quency and thus responds to the input at that frequency regard- 
less of the error of that input. 

This problem first manifested itself in the analysis of 
the Titan data discussed in the following section. A parallel 
effort with the three mass check case shed some insight on the 
effects of truncation. For example, two separate force time 
histories were applied to the model depicted in Figure 3. One 
force time history consisted of a decaying sine wave shown in 
Figure 4. The same decaying sine function with a superimposed 
steady-state value is shown in Figure 5. These two forcing 
functions characterize one of the most common truncation 
effects (Reference 4) . Figure 4 could represent a lateral 
force at a payload/booster interface extracted from recorded 
flight data. Figure 5 could be obtained in a like manner and 
represent a longitudinal interface load, i.e., in the direction 
of thrust. 

The forcing function shown in Figure 4 was used to drive 
the 3 mass check model and obtain the acceleration as a function 
of time for the middle mass at degree of freedom number 2 (see 
Figure 3). In addition, the Fourier transform of the forcing 
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Figure 5. Decaying Sine Forcing Function with 
Steady State 
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function in Figure 4 was generated using an FFT after the time 
history was "windowed" by a Tukey window (Reference 4) . Note 
that the time history in Figure 5 shows the results of the 
Tukey window application. This Fourier transform was then used 
to drive the model using equation (28) to calculate the response 
of degree of freedom number 2. The spectral results of the 
impedance analysis is shown in Figure 6. But for the true 
comparison this spectral response was inverse Fourier trans- 
formed to obtain a time history to correlate with the time 
domain response. This correlation is shown in Figure 7. 

The effect of the window is shown perhaps more in driving 
the model in the frequency domain with the force in Figure 5. 
Figures 8 and 9 show the Fourier spectral response of degree 
of freedom 2 with and without a Tukey window. Figures 10 and 
11 again depict a comparison of both sets of frequency domain 
response to the time domain response. It is obvious from these 
plots that for this case the Tukey window improves the response 
correlations . 

It is apparent from the results obtained that the low 
frequency environments can be predicted with the impedance 
technique. Better correlation with the time domain solutions 
could possibly be made with a different window or increasing 
the frequency resolution, however, enough correlation is 
apparent to attempt to use the technique for more complicated 
systems . 
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Figure 7. Time Domain Comparisons of Impedance 
Technique Responses and Time Domain 
Responses to Decaying Sine with no 
Steady State. 
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Figure 8. Fourier Spectrum of Response to Decaying 

Sine with Steady State but no Tukey Window. 
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Figure 9. 


Fourier Spectrum of Response tc Decaying 
Sine with Steady State with a Tukey Window. 
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Figure 11. Time Domain Comparisons of Impedance 
And Time Domain Solutions With A 
Tukey Window. 
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Analysis of the Titan Data 

The true test of the impedance technique was to analyze 
real available flight data. The objectives of this evaluation 
included discovering any pitfalls in the practicality of the 
overall impedance approach and to identify key items necessary 
for a payload user's guide. 

Available telemetered flight data from the first four 
Titan Centaur Launch Vehicle flights were digitized for use 
in this evaluation. Briefly, these flights can be described 
as follows: 

1. E-l generally consisted of a two stage Titan, solid 
rocket strap-ons. Centaur upper stage and Dynamic 
Simulator of the Viking payloads, with shroud; 

2. E-2 , same as E-l with a Helios payload; and 

3. E-3 and E-4 were the same as described above with 
Vikings A and B as payloads. 

The necessary information for solving the impedance equations 
(28) and (29) included the external forces for driving the 
analytical models and interface accelerations for comparison 
with the results. 

All recorded flight data, including the pulse code 
modulated (PCM) diagnostic data for the boost-vehicle as well 
as the FM/FM payload data were digitized at four hundred 
samples per second, and furnished by Langley Research Center. 

A further description of these flights and telemetered data 
can be found in Reference 5. 

Not all of the digitized data channels were used in this 
analysis. The transducers that were used measured chamber 
pressures for Stage I of Titan and acceleration (accelerometer 
number CM101A) near the Centaur /payload interface in the 
longitudinal axes (see Figure 12) . These measurements were 




representative of the applied external forces and interface 
accelerations for some of the transient events. In that regard, 
the event chosen for the impedance technique demonstration was 
the Titan Stage I burnout. This event by definition had 
fewer externally applied forces than say the liftoff event 
and it was past the major aerodynamic environments. Thus, 
the assumption that the engine pressure time histories 
describe sufficiently the external force/torque vector seems 
good for this event. 

As discussed in the previous section, some data manipula- 
tion is necessary for the flight data in order to obtain 
meaningful spectra. However, the windows discussed were not 
applied to the Titan data (they were developed later) although 
a similar technique was used. Linear ramps were multiplied 
with the data for fixed increments at the beginning and end 
of each extracted time history. The result was a reduction 
in the truncation errors much like using a "filter" type 
window. Figures 13 through 24 show the resulting time 
histories for engine pressure (denoted by TP sensor numbers) 
and acceleration for CM101A. The length of the time histories 
was 5.12 sec which corresponds to 2048 samples at 400 samples 
per sec. It can be seen in reviewing, say the pressure 
transducer TP3015, that slight differences exist from one 
flight versus another in not only amplitude but signature as 
well, a point which will be discussed in more detail later. 
Fourier spectra of these signals were calculated using an 
FFT routine and are shown in Figures 25 through 48. Notice 
each spectrum is plotted twice, the first showing the spectrum 
up to the Nyquist frequency 200 Hz, and the second one on one 
cycle log-log for more visual resolution in the lower frequen- 
cies. Here again, reviewing the spectra of any one transducer 
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Figure 13, 


Pressure Time History for TP3015 From 

E “1 Stage I Burnout With Front and End Ramps 
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Figure 15. 
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Acceleration Time History for CM101A From 
E-l Stage I Burnout With Front and End Ramps 
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Figure 19. Pressure Time History for TP3015 From E-3 
Stage I Burnout With Front and End Ramps 
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Figure 20. Pressure Time History for TP3016 From E-3 
Stage I Burnout With Front and End Ramps 
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Figure 21. Acceleration Time History for CM101A From 

E-3 Stage I Burnout With Front and End Ramps 
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Figure 22. Pressure Time History for TP3015 From E-4 
Stage I Burnout With Front and End Ramps 
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Figure 23. Pressure Time History for TP3016 From E-4 
Stage I Burnout With Front and End Ramps 
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Figure 24. Acceleration Time History for CM101A From 

E-4 Stage I Burnout With Front and End Ramps 
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Figure 25. Fourier Spectrum of TP3015 From E-l Stage I 
Burnout 
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Fourier Spectrum (lb-sec) (x.2248 N-sec) 



Figure 26, Fourier Spectrum of TP3015 From E-l Stage I 
Burnout (10-100 Hz) 
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Figure 27. Fourier Spectrum of TP3016 From E-l Stage I 
Burnout 
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Figure 28. 


Fourier Spectrum of TP3016 From E-l Stage I 
Burnout (10-100 Hz) 
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Figure 30. Fourier Spectrum of CM101A From E-l Stage I 
Burnout (10-100 Hz) 
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Fourier Spectrum (lb-sec) (x.2248 N-sec) 



Figure 31. Fourier Spectrum of TP3015 From E-2 Stage I 
Burnout 
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Fourier Spectrum (lb-sec) (x.2248 N-sec) 



Figure 32. 


Fourier Spectrum of TP3015 From E-2 Stage I 
Burnout (10-100 Hz) 
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Fourier Spectrum (lb-sec) (x.2248 N-sec) 



Figure 33. Fourier Spectrum of TP3016 From E-2 Stage I 
Burnout 
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Figure 34. 


Fourier Spectrum of TP3016 From E-2 Stage I 
Burnout (10-100 Hz) 
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Figure 35. Fourier Spectrum of CM101A From E-2 Stage I 
Burnout 
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Figure 36. Fourier Spectrum of CM101A From E-2 Stage I 
Burnout (10-100 Hz) 
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Figure 37. Fourier Spectrum of TP3015 From E-3 Stage I 
Burnout 
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Figure 38. Fourier Spectrum of TP3015 From E-3 Stage I 
Burnout (10-100 Hz) 
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Figure 39. Fourier Spectrum of TP3016 From E-3 Stage I 
Burnout 


60 


Fourier Spectrum (lb-sec) (x.2248 N- 


c 


a 

<u 

CO 


t 


c 



Figure 40 


Fourier Spectrum of TP3016 From E-3 Stage I 
Burnout (10-100 Hz) 
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Figure 41. Fourier . Spectrum of CM101A From E-3 
Stage I Burnout 
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Fourier Spectrum (g-sec) 
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Figure 42. Fourier Spectrum of CM101A From E-3 
Stage I Burnout (10-100 Hz) 
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Figure 43. Fourier Spectrum of TP3015 From E-4 
Stage I Burnout 
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Fourier Spectrun (lb-sec) (x.2248 N-sec) 



Figure 44. Fourier Spectrum of TP3015 From E-4 
Stage I Burnout (10-100 Hz) 
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Fourier Spectrum (lb-sec) ( x# 2248 N-sec) 



Figure 45. Fourier Spectrum of TP3016 From E-4 
Stage I Burnout 
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Fourier Spectrum (lb-sec) (x.2248 N-sec) 



Figure 46. Fourier Spectrum of TP3016 From E-4 
Stage I Burnout (10-100 Hz) 
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Figure 47. Fourier Spectrum of CM101A From E-4 
Stage I Burnout 
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Fourier Spectrum (g-sec) 
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Figure 48, Fourier Spectrum of CM101A From E-4 
Stage I Burnout (10-100 Hz) 
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from one flight to the next, the differences are quite 
evident . 

In order to generate the interface accelerations using 
the impedance equations, the analytical models for Stage I 
burnout of all four boosters and the payloads were resurrected. 
Booster model information was used for the complex transfer 
and point admittance coefficient matrices and the payload 
data generated the necessary complex payload impedance 
matrices . 

Since the structural configurations were the same for 
the boosters (with the exception of slight differences in 
propellant line dynamics), the only difference in model para- 
meters from one flight vs another was related to structural 
weight (see Tables 1 through 3) ; the largest difference 
being in the residual Stage I slosh fuel at burnout. Fifty 
modes were used for each booster in the analysis set. All 
modes above 50 Hz were truncated since the frequency regime 
of interest was approximately 18 Hz, the frequency content of 
a typical Stage I burnout. 

Tables 4 through 7 show a comparison of the modal 
frequencies for all longitudinal, pitch, roll and yaw modes. 

It is evident that in comparison the models vary only slightly 
in modal frequencies. The same is true for the eigenvectors 
used for the transfer characteristics. 

The modal damping used for the booster model was per- 
turbated a number of times to obtain what was felt to be 
optimum agreement with the first data spectra. Originally 
the booster damping chosen corresponded to that shown in 
Figure 49 which has been used in Titan loads analysis for 
some time. However, as stated before, the Stage I burnout is 
basically an 18 Hz event. The rise in the damping curve at 



TABLE 1. LONGITUDINAL WEIGHT SUMMARY FOR 
BOOSTER STAGE I BURNOUT MODELS 


Item 

E-l 

Wt. (N) 

E-2 

Wt. (N) 

E-3 

Wt. (N) 

E-4 

Wt. (N) 

Centaur 

1.580x10 5 

1.57 6x10 5 

1.576xl0 5 

1.570x10 

Interstage 

Adapter 

5.004xl0 3 

5.004xl0 3 

5.058xl0 3 

5.058x10 

Fairing 

3.205xl0 4 

3.205xl0 4 

3.205xl0 4 

3.205x10' 

Core (Struct) 

9.963xl0 4 

9.963xl0 4 

9. 963xl0 4 

9.963x10' 

Stage II Ox 

1.900x10 5 

1.921xl0 5 

1.913x10 5 

1.910x10 

Stage II Fuel 

1.072xl0 5 

1.077xl0 5 

1.0 64x10 5 

1.054x10 

Stage I Ox 

5.738xl0 2 

5.293xl0 2 

5.293xl0 2 

5.738x10 

Stage I Fuel 

1.066xl0 4 

1.108xl0 4 

2.584xl0 3 

3.679x10 

Total 

6.032xl0 5 

6. 10 2x10 5 

5. 9 50x10 5 

5.944x10 
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TABLE 2. LATERAL WEIGHT SUMMARY FOR BOOSTER 
STAGE I BURNOUT MODELS 


Item 

E-l 

Wt. (N) 

E-2 

Wt. (N) 

E-3 

Wt. (N) 

E-4 

Wt. (N) 

Centaur 

1.585xl0 5 

1.57 8x10 3 

1.578x10 5 

1.574x10 

Interstage 

Adapter 

5.004xl0 3 

5.004xl0 3 

5. 058x10 3 

5.058x10 

Fairing 

3.205xl0 4 

3.205xl0 4 

3.205xl0 4 

3.205x10 

Core (Struct) 

9.962xl0 4 

9.962xl0 4 

9. 962xl0 4 

9.962x10 

Stage II Ox 

1.900x10 5 

1.92 1x10 5 

1. 912xl0 5 

1.910x10 

Stage II Fuel 

1.072xl0 5 

1.077xl0 5 

1.0 64x10 3 

1.054x10 

Stage I Ox 

5.738xl0 2 

5.738xl0 2 

5.738xl0 2 

5.738x10 

Stage I Fuel 

1.066xl0 4 

1.107xl0 4 

2. 580x10 3 

3.679x10 

Total 

6.036xl0 5 

6.060xl0 5 

5.953xl0 5 

5.948x10 
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TABLE 3. SLOSH WEIGHT SUMMARY FOR BOOSTER STAGE I 

BURNOUT MODELS 


E-l 

Item Slosh Wt. (N) 

Centaur 

4 

LO 2 1.432x10 

LH 2 2. 153x10 3 

Stage II 

Ox 3. 674x10 4 

Fuel 2.285xl0 4 


E-2 

Slosh Wt. Oft 


1.432xl0 4 

2.153xl0 3 


3.674xl0 4 

2.285xl0 4 


E-3 

Slosh Wt. (N) 


1.432xl0 4 

2.153x103 


3.674xl0 4 

2.285xl0 4 


E-4 

Slosh Wt. (N) 


1.432xl0 4 
2 . 153xl0 3 


3.674x10 

2.285x10 


Stage I 


Ox 

8.896x10° 

8.896x10° 

8.896x10° 

8.896x10' 

Fuel 

7.967xl0 3 

8.380xl0 3 

1 . 957xl0 3 

2.576x10 
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TABLE 4. LONGITUDINAL MODAL FREQUENCIES FOR STAGE I 
BURNOUT BOOSTER MODELS 


Mode No. 

E-l 

Freq. (Hz) 

E-2 

Freq. (Hz) 

E-3 

Freq. (Hz) 

E-4 
Freq. (I 

1 

0.0 

0.0 

0.0 

0.0 

2 

18.1 

18.1 

18.4 

18.4 

3 

21.9 

21.8 

22.4 

22.3 

4 

26.8 

26.7 

27.2 

27.2 

5 

31.5 

31.5 

31.6 

31.5 

6 

37.4 

37.4 

37.4 

37.4 

7 

44.0 

44.0 

44.0 

44.1 
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TABLE 5. ROLL MODAL FREQUENCIES FOR STAGE I BURNOUT 

BOOSTER MODELS 


js 

[{■ Mode No . 

E-l 

Frea. (Hz) 

E-2 

Frea. (Hz) 

E-3 

Freq. (Hz) 

E-4 

Frea. (Hz) 

A : 

; 1 

0.0 

0.0 

0.0 

0.0 

2 

13.3 

13.3 

13.3 

13.3 

3 

15.9 

15.9 

15.9 

15.9 

4 

34.4 

34.4 

34.4 

34.4 

5 

38.1 

38.1 

38.1 

38.1 
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TABLE 

6. PITCH MODAL FREQUENCIES 
BOOSTER MODELS 

FOR STAGE I 

BURNOUT 

Mode No. 

E-l 

Frea. (Hz) 

E-2 

Frea. (Hz) 

E-3 

Freq. (Hz) 

E-4 

Freq. (Hz) 

1 

0.0 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

0.0 

3 

0.7 

0.7 

0.8 

0.9 

4 

0.9 

0.9 

0.9 

0.9 

5 

1.0 

1.0 

1.0 

1.0 

6 

1.0 

1.0 

1.0 

1.0 

7 

1.2 

1.2 

1.2 

1.2 

8 

3.1 

3.1 

3.1 

3.1 

9 

6.5 

6.5 

6.5 

6.5 

10 

10.5 

10.5 

10.4 

10.5 

11 

13.4 

13.4 

13.4 

13.4 

12 

16.0 

16.0 

16.0 

16.0 

13 

18.2 

18.1 

18.2 

18.2 

14 

28.7 

28.7 

28.7 

28.7 

15 

34.9 

34.9 

35.0 

35.0 

16 

36.7 

36.7 

36.8 

36.8 

17 

40.0 

40.0 

40.1 

40.1 

18 

46.2 

46.5 

46.2 

46.2 

19 

48.1 

48.0 

48.0 

48.0 
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TABLE 7 . YAW MODAL FREQUENCIES FOR STAGE I BURNOUT 

BOOSTER MODELS 


Mode No, 

E-l 

Frea. (Hz) 

E-2 

Freq. (Hz) 

E-3 

Freq. (Hz) 

E-4 

Freq. (Hz) 

1 

0.0 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

0.0 

3 

0.7 

0.7 

0.8 

0.9 

4 

0.9 

0.9 

0.9 

0.9 

5 

1.0 

1.0 

1.0 

1.0 

6 

1.0 

1.0 

1.0 

1.0 

7 

1.2 

1.2 

1.2 

1.2 

8 

3.1 

3.1 

3.1 

3.1 

9 

6.5 

6.5 

6.5 

6.5 

10 

10.5 

10.5 

10.5 

10.5 

11 

13.4 

13.4 

13.4 

13.4 

12 

16.0 

16.0 

16.0 

16.0 

13 

18.2 

18.2 

18.2 

18.2 

14 

28.7 

28.7 

28.7 

28.7 

15 

34.9 

34.9 

35.0 

35.0 

16 

36.9 

36.9 

37.0 

37.0 

17 

40.0 

40.0 

40.1 

40.1 

18 

46.2 

46.5 

46.2 

46.2 

19 

48.1 

48.0 

48.0 

48.0 
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approximately this frequency appeared to attenuate the response 
to an extreme degree. Therefore, the optimum modal damping for 
good comparison to flight data spectra of measurement CM101A 
was chosen as 1.5% for all modes. Since there appears to be 
very little justification for using the curve in Figure 49 
aside from precedent, this assumption seems valid. 

The payload models varied understandably from one flight 
to another with the exception of the E-3 and E-4 flights which 
were, of course, both Viking payloads. A detailed breakdown of 
each payload’s dynamic characteristics would be beyond the 
scope of the intent of this report. Instead, frequency domain 
plots of the impedances for each payload in the longitudinal 
direction are shown in Figures 50 through 53. Table 8 gives 
total weights and centers of gravity for comparisons. The 
impedance shown are based on mode shapes truncated above 
50 Hz and include the individual payload damping. For the 
Viking models, the damping used was based on coupled damping 
from modal coupling discarding the off-diagonal terms. 

Helios assumed a damping curve as a function of frequency 
identical to that shown previously in Figure 49, and the VDS 
damping was obtained directly from the modal survey of the 
flight structure. 

The Titan data analysis was broken into two parts. One 
part of the effort demonstrated the solution for the payload 
interface accelerations using a computer program written to 
solve equation (28). For each of the four flights, the 
Fourier spectra of the engine pressure time histories were 
used to drive the models resulting in the interface accelera- 
tion spectra in the longitudinal direction. These calculated 
Fourier response spectra are shown in Figures 54 through 57. 
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Figure 50. Longitudinal Payload Impedance For VDS. 
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Impedance (lb/g) (x.2248 N/g) 



Figure 51. 


Longitudinal Payload Impedance For Helios. 
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TABLE 8. COMPARISON OF PAYLOAD TOTAL WEIGHTS 
AND C.G. LOCATIONS 


Payload 

Total 
Weight (N) 

X M 

C.G. Location* 

Y CM) Z CM) 

VDS 

3.48xl0 4 

3.64x10° 

3. 35x10“ 2 -1.52x10" 2 

Helios 

1.56xl0 4 

2.13x10° 

-3.95x10” 3 -3.13xl0" 4 

Viking A 

3.42xl0 4 

3.33x10° 

2.08xl0” 2 -1.83xl0” 2 

Viking B 

3.42xl0 4 

3.33x10° 

2.08x10” 2 -1.83xl0” 2 


* C.G. locations are referenced to the booster interface with 
a coordinate system shown below: 
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Figure 54. Comparison of CM101A Flight Data 

Spectrum With Impedance Analytical 
Predictions for E-l Stage I Burnout 
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Figure 55. Comparison of CM101A Flight Data 

Spectrum With Impedance Analytical 
Predictions for E-2 Stage I Burnout 
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Figure 56. Comparison of CM101A Flight Data 

Spectrum With Impedance Analytical 
Prediction for E-3 Stage I Burnout 


87 


Fourier Spectrum (g-sec) 





The second part of the analysis involved using the flight 
data to solve equation (29). The Helios flight, E-2 , was 
chosen as the basis for the original interface data from which 
the low frequency environments were predicted. Most of the 
modes of the Helios payload can be considered to be longi- 
tudinal in nature based on information from previous work. 

For this reason, most of the response shown in the CM101A 
measurement during the E-2 flight should be due to longitudinal 
motion and very little to lateral and bending mode feedback. 

A separate computer program was written to solve equation (29) 
and the results for the predicted Viking and VDS longitudinal 
environments are shown in Figures 58 and 59. 

As expected the first obvious problem encountered in 
using the CM101A data pertained to its location. The models 
used for the generation of the booster admittances, i.e., in 
this case the transfer characteristics of the coupled Centaur 
and two stages of the Titan III, and the structural response 
measured by CM101A were inconsistent. CM101A was physically 
located on a forward equipment truss at the outer circum- 
ference of the Centaur structure. In that location, the 
accelerometer picked up not only the desired longitudinal 
vehicle accelerations but any bending mode responses as well 
as local secondary structure responses. 

One other model/flight data correlation problem 
encountered was due to the Titan second stage ignition. 

Close comparisons of engine thrust time histories and 
accelerometer time histories in Figures 13 through 24 show 
a response on the accelerometer that does not appear on the 
pressure time histories. This response occurs after the first 
stage shutdown, i.e., where the acceleration bagins to pick up 
a dc level, and contains 20 to 22 Hz frequency oscillations. 
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Figure 58. Comparison of CM101A E-l Spectrum With 
Impedance Ratio Response From E-2 
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Fourier Spectrum (g-sec) 



Figure 59. Comparison of CM101A E-3 Spectrum With 
Impedance Ratio Response From E— 2. 
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The Stage I burnout models, however, have neither the true 
modeling of the physical separation of the first and second 
stages nor the forcing function corresponding to the second 
stage fire. For these reasons close comparisons of flight 
data and analytical predictions are not expected above the 
18 Hz first stage burnout. It should also be pointed out 
that the second stage ignition problem will be different from 
one flight to another. For this reason in particular, the 
frequency content of the spectra will be different also. 

However, even with these discrepancies, the correlations 
of flight data and analytical prediction are quite good 
(Figures 54 through 57) . For the desired mode of response 
(approximately 18 Hz) , the technique appears to correlate in 
amplitude and bandwidth. With proper signal conditioning and 
modeling the correlation at the other frequencies can be 
improved. These points are important for consideration of a 
payload user’s guide. 



Application to STS Payload 

As a further demonstration of the impedance technique the 
second phase of the analytical effort was directed at predicting 
the low frequency environment for an existing STS payload. The 
primary objective of this analysis was to compare results with 
other current environmental predictions from state-of-the-art 
techniques. It is hoped that the results of these predictions 
can then be used for comparisons with flight measured data 
from the development flights. 

The models used for this set of analyses were chosen to be 
consistent with those used for the loads cycles done by the 
STS project. They were obtained through LRC and were 
comprised of eigenvalues and eigenvectors for the STS 
finite element model of the liftoff and landing configurations. 
Two hundred eigenvectors were used for each configuration 
which included the models of two of the three payloads in the 
payload bay for this flight, the DFI and SP . All modes for 
each loads case were assumed to have a modal damping ratio of 1%. 

The third payload in the bay, the Long Duration Exposure 
Facility (LDEF) , was chosen as the payload for these analyses. 
Forty eigenvalues and eigenvectors were again furnished by LRC. 
The discrete mass matrix for the LDEF, necessary for the genera- 
tion of the inertial coupling coefficients matrix on the right- 
hand side of equation (22) was not available and had to be 
generated from discrete model data with in-house finite 
element routines (see Figure 60) . Orthogonality checks were 
calculated to ensure that the modes and regenerated mass 
matrix were consistent. Like the booster models, the LDEF modes 
were assumed to have 1% damping throughout. 

The agreed upon flight events for analysis were liftoff 
and landing. The liftoff event analyzed was defined as a 
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"high performance fast timing" engine thrust with a mismatch 
between right and left SRBs . The SRB mismatch in this case 
was 0.098 sec. 

A total of 62 discrete forcing functions were used to 
drive the analytical model for liftoff. Each defined time 
history was 10.0 sec in length. To obtain the Fourier 
spectrum of each force, they were first digitized, i.e., 
linearly interpolated, at 102.4 samples per second. This 
resulted in a total number of samples of 1024 with a Nyquist 
of approximately 50 Hz. The same Tukey window discussed in 
the previous sections was used for the liftoff forcing func- 
tions. Representative Fourier spectra of these forces were 
plotted and are shown in Appendix A. 

The spectral results of the impedance predictions for the 
landing event are shown in Figures 61 through 66. These plots 
show the frequency bandwidths of maximum and minimum interface 
acceleration for the LDEF during the LP501R liftoff event. 

For a check with previous predictions, the complex Fourier 
spectra were inverse transformed to obtain a time history for 
amplitude and frequency content comparisons. The inverse 
Fourier transform time histories are presented in Figures 67 
through 72. Shown also on the plots are the corresponding 
peak accelerations resulting from a time domain analysis. 
Corresponding, Z, or pitch, responses for the aft end of the 
LDEF were not available for comparison. 

The Tukey window effect on the resulting time histories 
is best shown in Figure 67 which is the longitudinal accelera- 
tion at a forward LDEF attach point. Because of the window, 
this response begins and ends at approximately zero amplitude. 
However, since this windowing technique of signal conditioning 
does not alter the signal in the remaining areas, its effect 
can be ignored. 
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Figure 61. Fourier Spectrum of X Response At 
Forward Attach Point Right Hand 
Side For Lift-Off. 


Fourier Spectrum (g-sec) 



Figure 62. Fourier Spectrum of Z Response At 
Forward Attach Point Right Hand 
Side For Lift-Off. 
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Figure 63. Fourier Spectrum of X Response At 
Forward Attach Point Left Hand 
Side For Lift-Off. 
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Figure 64. Fourier Spectrum of Z Response At 
Forward Attach Point Left Hand 
Side For Lift-Off. 
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Figure 65. Fourier Spectrum of KEEL Response 
For Lift-Off. 
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Figure 66. Fourier Spectrum of Aft Z Response 
For Lift-Off. 
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Figure 68. Response Z Forward Attach Right 
Hand Side For Lift-Off. 
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Figure 69* Response X Forward Attach Left 
Hand Side For Lift-Off* 
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Figure 70. Response Z Forward Attach Left 
Hand Side For Lift-Oft, 
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Figure 72. Response Z Aft For Lift-Off. 
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The correlation of the impedance technique results and 
the time domain solutions show generally good comparisons. 
Transcribing the peaks of the time domain accelerations to the 
inverse transform plots is agreed not to give conclusive 
comparison parameters. However, the peaks do appear to occur 
at the same points in time for both cases, indicating a minimal 
phase error. The largest discrepancy is shown in amplitude 
comparisons. Even in this area some of the amplitudes cal- 
culated with the impedance technique compare quite well. The 
best comparisons are shown in the longitudinal accelerations 
(Figures 67 and 69) . 

The final event analyzed with the impedance technique was 
a landing case. For this event 266 Fourier spectra were cal- 
culated for discrete forcing functions for the orbiter (repre- 
sentative loads and spectra can be found in Appendix B) . The 
event can generally be described as a high angle of attack with 
symmetric landing conditions. 

The problems encountered in analyzing the landing event 
resulted in defining analytical requirements for the use of the 
impedance technique for low frequency environment predictions. 
The analysis of this event pointed out the pitfalls of fre- 
quency resolution. The time histories of the forcing function 
for landing were of two different lengths. These are two 
basic external forces on the orbiter for landing (Appendix B) : 
the landing gear strut loads; and the aerodynamic loads due to 
such things as ground effects, etc. The strut force time 
histories were 2.0 seconds in length and all other loads were 
only 0.8 seconds in length. Since the behavior of the aero- 
dynamic loads were not known past 0.8 second and since it 
appeared that all oscillatory characteristics of the strut 
forces had diminished by 0.8 sec, the decision was made to 
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truncate the strut forces at 0.8 seconds for the analysis. 

The resulting sample rate based on the 512 data points used was 
640 samples per second with a Nyquist frequency of 320 Hz. 

The important thing of note, however, was that the fre- 
quency resolution for the 0.8 second signal was only 1.25 Hz. 
This resolution would be sufficient for describing character- 
istics of higher frequencies. However, the orbiter model 
contains low frequency bending modes at or around 5 Hz. With 
the assumed 1% modal damping ratio, this mode may have a band- 
width of 0.1 Hz. Thus, it is obvious that frequency descrip- 
tions of the input forces with a frequency resolution of 
1.25 Hz could and did give erroneous answers. 

The decision was made to discard the aerodynamic loads 
and perform the analysis with the full 2.0 second time history 
of the strut loads only. This decision of neglecting the 
aerodynamic forces was based on the general assumption that 
the aerodynamic forces contribute only a small part to the 
low frequency environment at the payload interface. The 
resulting resolution with only the strut forces becomes 
0.5 Hz, which comes somewhat closer to the desired 0.1 Hz. 

The results of the analysis with only the strut loads are 
shown in Figures 73 through 84. Here again the peaks from the 
time domain solution were transcribed to the plots for com- 
parison. 

These comparisons were not expected to be as good as the 
liftoff case generally because of the absence of aerodynamic 
loads. The worst comparisons are more evident in the "Z" 
direction acceleration, as expected, since this is the 
direction of the general aerodynamic forces. Also, the correla- 
tion seems to get worse from about 0.3 to 0.4 seconds on. This 
■again is when the aerodynamic forces are the most active. 
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Figure 73. Fourier Spectrum of X Response At 
Forward Attach Point Right Hand 
Side For Landing Due To Strut 
Forces Only. 
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Figure 74. Fourier Spectrum of Z Response At 
Forward Attach Point Right Hand 
Side for Landing Due to Strut 
Forces Only . 
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Figure 75. Fourier Spectrum of X Response At 
Forward Attach Point Left Hand 
Side For Landing Due To Strut 
Forces Only. 
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Figure 76. Fourier Spectrum of Z Response At 
Forward Attach Point Left Hand 
Side For Landing Due To Strut 
Forces Only. 
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Figure 77. Fourier Spectrum of KEEL Response For 
Landing Due to Strut Forces Only. 
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Figure 78. Fourier Spectrum of Aft Z Response For 
Landing Due to Strut Forces Only. 
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Figure 80. Response Z Forward Attach Right Hand 
Side For Landing Due To Strut Forces 
Only. 
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Figure 83 . Response At KEEL For Landing Due To 
Strut Forces Only. 
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Figure 84. Response Z Aft For Landing Due To 
Strut Forces Only. 
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DISCUSSION OF RESULTS 


The results of analyzing the Titan III-E data made some 
significant contributions to the evaluation of the impedance 
technique. The most important result was somewhat obvious from 
the beginning of the study. Extreme care has to be taken when 
recording and analyzing flight data for use in load prediction 
techniques. Transducers used for this phase of the study were 
not optimum for good predictions of the low frequency environ- 
ments or for description of the vehicle transfer characteristics. 

Besides pointing out data recording requirements, the 
Titan data analysis results consisted of new efforts in signal 
conditioning and spectral analysis. Efforts using ramps, 
windows, and FFTs refined the overall philosophy of the 
impedance technique. Even with these conditioning techniques, 
the discrepancies and extraneous signals on the flight data 
make evaluation of the spectral responses against the flight 
data quite difficult. 

These results do, however, point out one important point: 
determining the low frequency environment from a pure spectral 
standpoint is quite difficult. As the results of the STS 
analysis show, additional information is gained from the history 
obtained from the inverse transforms. 

The analysis of the liftoff event for the LDEF by far 
showed the best results. Even these results can be improved, 
however. As discussed in the previous sections, the differences 
in coupled damping versus the resulting damping of the impedance 
technique will never go away. The differences in the descrip- 
tion of the equations of motion at the interface that depend on 
the boundary conditions will always cause that result. However, 
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with in-flight information of the true transfer characteristics, 
these problems could be minimized. 

The set of analyses performed with the impedance technique 
for the liftoff case was done somewhat incorrectly to evaluate 
the impact on the results. The steady-state winds and gust 
loads (Appendix C) are represented by a dc level with a ,T l-cos" 
gust function. In normal time domain solutions the necessary 
initial conditions are derived to place the vehicle in equilib- 
rium before the time domain solutions begin. In the spectral 
analysis of the impedance technique, these forces were applied at 
time zero, forcing the model to respond to an "imaginary" tran- 
sient at time zero. Therefore, the responses shown will be changed 
when the winds are treated correctly. 

The landing event analysis probably had the most meaningful 
results for the intent of this study. These results when 
compared to the time domain analysis predictions show not only 
the effects of the aerodynamic forces but more importantly 
give insight into the errors caused by inadequate frequency 
resolution. 

Neither the liftoff or landing results shown should be 
used in any manner for design purposes for other payloads. The 
models used in these analyses are not necessarily the latest 
configurations. In addition, assumptions made in the impedance 
technique analysis, e.g., neglecting aerodynamic loads, cause 
the results to be somewhat unconservative. 

Based on the results of all of the analyses performed, a 
limited evaluation of the technique for the use of future 
payloads can be made. It appears that this technique can be 
used for any coupled payload system in a very inexpensive 
manner. The ratioing of the acceleration from one flight to 
the next, as developed in the equation of motion and demonstrated 
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with the flight data from Titan, probably has the least use for 
STS payloads due to the variety of interface locations in the 
bay. It may be valuable, however, in evaluating design changes 
once one loads analysis is conducted. 

The best use can be to achieve the removal of an integra- 
tion role with a final set of eigensolutions . This can be 
achieved in much the same way as was done in this study. In 
the future, payload 1 projects can obtain a data tape from the 
STS project that contains the best STS model and flight event 
cases. The payload designer can then perform the analysis to 
evaluate the design. This philosophy merely requires continual 
update of the STS dynamic characteristics and the latest flight 
data information. 

Having calculated the predicted low frequency environment 
for a particular payload, the evaluation can be made for the 
necessity of a more extensive loads and/or response analysis 
for the payload based on the spectral content of the interface 
environment and the payload impedance. If there is concern 
about the design margin, the impedance technique programs should 
give some insight into the modal degrees of freedom of the 
generalized forces that are contributing most to the environ- 
ment at that frequency bandwidth. This information should 
then be used to reduce the size of the models involved and 
again cut costs. 

In summary, based on the results of this study, with 
reusable boosters and an environmental data bank, payload 
organizations can use the impedance technique to: 

1. eliminate the necessity of integrated coupled 
analysis ; 

perform their own low cost environment predictions; 


2 . 



3. reduce analytical effort with spectral evaluation of 
coupled response; and 

4- in the event more detailed analysis is necessary, 
reduce the size of the dynamic models. 
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CONCLUSIONS 


An analytical technique has been derived that can be used 
to predict low frequency environments for payloads and at the 
same time lower the cost of the necessary analytical effort to 
calculate those predictions. Results of the analytical efforts 
during this study show that the frequency domain impedance 
technique can be used to define the interface accelerations, 
in both spectral and time history form, with sufficient 
accuracy to evaluate the payload design. By eliminating the 
necessity for final coupled eigensolutions , performing the 
response analysis in the frequency domain, and using the 
FFT to obtain Fourier spectra, the impedance technique reduces 
not only analytical integration effort, but computation costs 
as well. 

Future payloads developers using common boosters will be 
able to perform their own environmental prediction analyses 
and therefore possibly impact the design earlier in the payload 
program. By evaluation of the low frequency environments 
versus the payloads dynamic characteristics and design margin, 
early decisions for more or less extensive analyses can be 
made . 

Martin Marietta Corporation 
Denver Division 
Denver, Colorado 80201 
June 26, 1978 
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Figure A-l. Typical Fourier Spectrum of "l-cos" 

Wind Load Used For Lift-Off Analysis. 
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Fourier Spectrum (lb-sec) (x.2248 N-sec) 
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Figure A-2. Typical Fourier Spectrum of Orbiter 

Engine Thrust Used For Lift-Off Analysis. 
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APPENDIX B 
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Figure B-2. Typical Fourier Spectrum of Strut Loads 
From 0.8 sec Time Histories Used For 
Landing Analysis. 


Fourier Spectrum (lb-sec) (x.2248 N-sec) 



Figure B-3. Typical Fourier Spectrum of Strut Loads 
From 2.0 sec Time Histories Used for 
landing Analysis. 
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